Self-similar cosmological solutions with a non-minimally coupled scalar field 
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We present self-similar cosmological solutions for a barotropic fluid plus scalar field with Brans- 
Dicke-type coupling to the spacetime curvature and an arbitrary power-law potential energy. We 
identify all the fixed points in the autonomous phase-plane, including a scaling solution where 
the fluid density scales with the scalar field's kinetic and potential energy. This is related by a 
conformal transformation to a scaling solution for a scalar fleld with exponential potential minimally 
coupled to the spacetime curvature, but non-minimally coupled to the barotropic fluid. Radiation 
is automatically decoupled from the scalar field, but energy transfer between the field and non- 
relativistic dark matter can lead to a change to an accelerated expansion at late times in the 
Einstein frame. The scalar field density can mimic a cosmological constant even for steep potentials 
in the strong coupling limit. 
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I. INTRODUCTION 

Scalar fields have come to play a dominant role in re- 
cent years in theoretical models of the universe. This has 
usually been in the context of inflationary models of the 
very early universe where the self interaction potential 
energy density V{4') remains undiluted by the cosmolog- 
ical expansion. If the potential is sufficiently flat this 
can lead to an effective cosmological constant which can 
drive an accelerated expansion. The detailed nature of 
the evolution is driven by the specific form of the scalar 
field's potential energy. 

More recently observational evidence has been grow- 
ing that the energy density of the universe today may be 
dominated by a homogeneous (smooth) component with 
negative pressure - christened "quintessence" |^] . A suit- 
ably small cosmological constant A provides a good fit to 
the observational data , but the required low energy 
scale, A ^ (lO^'^eV)'*, seems unnatural in any supergrav- 
ity model. This has lead to an explosion of interest in the 
late-time evolution of self-interacting scalar fields, where 
radiation and matter contribute the dominant energy 
density at earlier times, but slowly rolling scalar fields 
may provide a dynamical mechanism to achieve a small 
effective cosmological constant today by mimicking a de- 
caying effective-A. The late time evolution of sufficiently 
steep potentials has been shown to be largely insensitive 
to initial conditions Q (though it should be remembered 
that it the case of a true cosmological constant there is 
also no dependence upon initial conditions!). This steep 
potential must become flat by the present in order to act 
as an effective cosmological constant today and there re- 
mains the "fine-tuning" problem of why the scalar field 
density comes to dominate at just the present epoch. 

Scalar fields with exponential potentials provide a 
strikingly simple model in which the energy density of 
the scalar field decreases at exactly the same rate as a 
barotropic fluid. Wetterich was the first to point out 
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that the energy density of a scalar fleld with an exponen- 
tial potential does not simply die away if the potential is 
too steep to drive inflation, but rather approaches an at- 
tractor solution Q where it mimics the equation of state 
of a barotropic fluid (radiation or dust) and remains a 
flxed fraction of the total energy density. We will refer 
to these solutions as self-similar or sca/m^l^ cosmological 
models where the dimensionless density parameter for 
each component, Q.i = Pi/pcrit, is a constant. 

Self-similar homogeneous cosmological solutions are in- 
variant under a global conformal rescaling of the metric. 
Such models may be expanding, but the physical state 
at different times differ only by a change in the over- 
all length scale In a Friedmann- Robertson- Walker 
(FRW) cosmological model this is equivalent to a rescal- 
ing of the cosmic time, 



t 



Tt 



(1.1) 



where F is a constant. Self-similarity requires that the 
evolution of the scale factor a{t) is a power-law, and the 
Hubble constant H = a/ a cx 1/t. 

We will consider the evolution of spatially flat FRW 
cosmologies containing a fluid, with density p, and pres- 
sure P, and a scalar fleld (p with self-interaction poten- 
tial V{(j)). In general relativity, the Friedmann constraint 
equation requires 



SttG 



N 



(1.2) 



where Gn is Newton's constant. Equation (IT) then re- 
quires that the matter content is invariant under a rescal- 



*Note that the term "scaling solution" has been applied else- 
where in the literature to describe any power-law evolution 
for the scalar field density in a radiation or matter dominated 
background. 
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ing p — > T~^p. Requiring P — > r~^P for all V then leads 
to a barotropic equation of state P = (7 — l)p where 7 
is a constant. Thus the familiar radiation (7 = 4/3) or 
pressureless matter (7 = 1) dominated FRW models can 
be described as self-similar solutions with 



(1.3) 



A scalar field has kinetic energy 0^/2 which only allows 



a constant shift of the scalar field, 



if we are 



to require (fp' ^ V '^(jp. In order to obtain a self-similar 
solution we require in addition 



y(0)^F(0 + 0r)=r-2y(0) 



(1.4) 



which is only compatible with a non-interacting field, 
V{(f>) = 0, or an exponential potential 



V{(j)) = Foexp {~Xh 



(1.5) 



where — SttGn, A is a dimensionless constant and we 
have K0r = (2/ A) InT. 

FRW solutions for scalar fields with exponential poten- 
tials where the kinetic energy and potential energy of the 
field remain proportional were proposed as a model for 
power-law inflation in the early universe by Luchinn and 
Matarrese and are the late-time attractor solutions 
(in the absence of other matter) for A^ < 6 ||l^,|ll[ . More 
recently attention has focused on the possible late-time 
evolution of scalar fields in FRW cosmologies contain- 
ing matter. Self-similar solutions are known for scalar 
fields with exponential potentials whose energy density 
scales with that of a barotrop ic fl uid yielding the same 
time dependence given in Eq. (1.3) for a barotropic fluid 
dominated solution ||^. These scaling solutions are the 
unique late-time attractors for sufficiently steep poten- 
tials (A^ > 37) 1^. Such a scalar field is so successful at 
scaling with the barotropic matter that the scalar field 
never comes to dominate the cosmological dynamics. A 
minimally coupled field can only mimic a cosmological 
constant at late times if the scale invariance of the scalar 
field is broken, leading to the deviation from the pure 
scaling solution which must be fine-tuned to occur at 
just the present epoch. 

In this paper we will consider the effect of non- minimal 
coupling of the scalar fleld to the spacetime curvature. In 
the next section we show that self-similar cosmological 
solutions are possible for scalar fields with simple power- 
law potentials if the scalar field has a Brans-Dicke-type 
coupling to the spacetime curvature. The equations of 
motion reduce to an autonomous phase-plane whose fixed 
points correspond to self-similar solutions. We show that 
this system is conformally related to a scalar field with ex- 
ponential potential and minimally coupled to the space- 
time curvature, but with an explicit coupling to matter, 
studied in Refs. p^-p^. This coupling is sensitive to the 
fluid equation of state and can produce a natural change 
in the nature of the solution after matter-radiation equal- 
ity. In contrast with the case of a minimally coupled 



scalar field there is the possibility of recovering an effec- 
tive late time cosmological constant even if the potential 
remains steep , depending on the value of the dimen- 
sionless non-minimal couplin g co nstant and we investi- 
gate this possibility in section 
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II. BRANS-DICKE-TYPE COSMOLOGY WITH 
POWER-LAW POTENTIAL 

It is possible to obtain self-similar solutions for scalar 
fields with arbitrary power-law potentials , V{(j>) — 
Vo(k^)^", if we go beyond Einstein's theory of general 
relativity and allow the field to be non-minimally coupled 
to the spacetime curvature. In Brans-Dicke gravity p^ ] 
with dimensionless parameter w, Newton's constant is 
replaced by a dynamical field G = w/(27r0^) and the 
generalised Friedmann equation requires 



Stt lu 
y 2^ 



^^^ + V{^)-^H^^] . (2.1) 




In t he o riginal Brans-Dicke theory where V{(l)) = 0, 
Eq. (1.1) is compatible with a global re-scaling^ of the 
barotropic fluid density p — s- Tpp and ((> — > T^(f) for all 
Tp/T'^ — r~^. The self-similar solution for a barotropic 
fluid in Brans-Dicke gravity was given by Nariai 
In the presence of a power-law potential, V((t>) 
Vb(K0)^", (absent in the pure Brans-Dicke theory [ |18[ ) 
we require in addition that V{4>)/4>^ T~''^\V {(f)) / (fP'] 
which specifies r~i/("-i)(/) and p r-2"/(»-i)p. 
Thus self-similar solutions can exist for "quintessence" - 
type power-law scalar field potentials [^,g3|, so long as 
this scalar field has Brans-Dicke type coupling to the 
spacetime curvature. 

Motivated by the preceding qualitative discussion of 
self-similar solutions we examine the cosmological evolu- 
tion of one of the simplest forms of non-minimal coupling 
proposed by Brans and Dicke |l^ which is described by 
a single dimensionless parameter oj. In addition to the 
original Brans-Dicke Lagrangi an pq | , we introduce a self 
interaction potential V{(f)) ||24|-|27|. The action is 

^ 9^ ^matter 

(2.2) 



where upper/lower signs should be chosen to ensure 
ibo; > and hence a positive gravitational coupling. The 
self-interaction potential is taken to be a power-law, 



^The global rescaling leaves the dimensionless Brans-Dicke 
parameter, lo, invariant. A local rescaling V{(f)) rescales the 
Brans-Dicke parameter |19|. 
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(2.3) 



This action can be re-expressed as a theory of interact- 
ing matter fields in general relativity with fixed Newton's 
constant Gn [Hi if we define quantities in the confor- 
mally related Einstein frame with respect to the rescaled 
metric 



(2.4) 



The scalar field is now minimally coupled to metric Qfiu^ 
but non-minimally coupled to the other matter fields. 
We will re-express the scalar field in terms of a field 
which has a canonical kinetic term in the Einstein frame, 
which requires that 



Kcp = exp 



Kip 



v/2(3 + 2uo) 



(2.5) 



The scalar field will have a non-negative energy density 
in the Einstein frame so long as uj > —3/2. Thus we 
consider to with any value > —3/2, which includes, for 
example, the coupling of the string-theory dilaton where 
= [|9|. 

The effective self-interaction potential for canonical 
field ip in the Einstein frame is of exponential form and 
given by (see, e.g., Refs. E^M) 



where Vq = (4u;)^Vb and 



F(0) = Vbexp(-AKV) , (2.6) 



A = 



3 + 2a; 



(2-n) 



(2.7) 



A. Autonomous Phase— Plane 

Using the Hubble rate, H, fluid density, p, and scalar 
field density, = {l/2){d'p/ dtf + V{ip), defined in the 
Einstein frame, the evolution equations are 



The non- minimal coupling of the scalar field in the Brans- 
Dicke frame leads to an explicit energy transfer in the 
Einstein frame, between the scalar field and the fluid. 



dQ 
di 



dip P-3P 



di ^2(3 + 2lu) 



(2.12) 



The same system of equations was considered, starting 
from a different motivation, by Wetterich [ p!3| , and more 
recently by Amendola jl4|. In the case dQ/dt = 0, these 
equations reduce to previous studies of minimally coupled 
scalar fields with exponential potentials 
Following Ref. M, we define. 



dp 
di 



(2.13) 



The evolution equations for a barotropic fluid where P — 
(7 — l)p can then be written as [pilj] 



x' - -3a; + ^y^ + ^x[2x^ + 7(1 - x^ ~ y^)] 

+ W{l~x^ ~y^), (2.14) 

y' = ->^\ll^y + \y[^^^ + 7(1 - 2;' - y')], (2.15) 

where a prime denotes the derivative with respect to = 
/ Hdi, the logarithm of the scale factor in the Einstein 
frame, and we have used the constraint equation. 



^+X^+y^ = l 



(2.16) 



We parameterise the energy transfer in terms of the 
dimensionless quantity W , where 



W 



3 + 2cj 



4-37 



(2.17) 



Equations (2.14) and (2.15) define an two-dimensional 



autonomous phase-plane whenever W can be written as 
function of x an d y. The Brans-Dicke-type theory de- 
fined in Eq. (2.2) naturally leads to an interaction with 
constant W , although the analysis could be extended to 
more general W{x^y) |l5|Jl6|l . 

The constraint Eq. (|2.16| ) restricts physical solutions 
with non-negative fluid density to < -I- j/^ < 1 
and so the evolution is completely described by trajecto- 
ries within the unit disc. In the following discussion we 
will only consider expanding cosmologies corresponding 
to the upper-half disc, y > 0. 



subject to the Friedmann constraint. 



B. Self-similar solutions 



3 r 2\dt 



V{iP) 



(2 11) Fixed points of the system, [xi^yP) where x' — Q and 

y' — 0, correspond to power-law solutions for the scale 



3 



factor and logarithmic evolution of the scalar field with 
respect to the cosmic time in the Einstein frame: 



d (X , Kip (X qlnt , 
where the constants p and q are given by, 

2 

P 



(2.18) 



q 



37(1 + 
37(l-a;f-y2) + 6a:f 



(2.19) 



These are self-similar solutions where the dimcnsionless 
density parameter for the barotropic fluid in the Einstein 
frame, given from the constraint Eq. (2.16) 



9 - 



1 



(2.20) 



is a constant and therefore so is the dimcnsionless density 
of the scalar field 



ri,/, — 



(2.21) 



The scalar field has an effective barotropic index given 

by, 



14: 



y 



(2.22) 



The existence and stability of critical points, already 
defined in the Einstein frame, remain unchanged under a 
conformal transformation back to the Brans-Dicke frame. 
The evolutionary behaviour for the scale factor and scalar 
field is however modified. If we denote the evolutionary 
behaviour of the scale factor a, and scalar field (j) in the 
Brans-Dicke frame by, 



i9 



(2.23) 



exponents can be related to those in the Einstein frame 

by, 



P = P + 



q = 



ip-i)q 



v/2(3 + 2cj) - g ' 
2q 



^y2{3 + 2uj)-q 



(2.24) 



Depending on the values of the parameters A, 7 and W 
we can have up to five fixed points in the Einstein frame. 
The nature and stability of each point in the phase-plane 
is briefly reviewed below. A full analysis of the existence 
and stability is given in the Appendix. See also Ref. [n4| . 



1. 2-Way, Matter-Kinetic scaling solution 
2W 



Xi 



3(2 - 7) 



yi = 0. 



(2.25) 



This solution lies on the a:-axis where the scalar field 
potential is negligible, and the scalar field's density in 
the Einstein frame is dominated by its kinetic energy, 
leading to a stiff equation of state for the scalar field, 
7^, = 2. We have a power-law solution of the form given 
in Eq. ( ^.18| ) with 



Pi 



qi 



6(2-7) 



97(2-7)-t-4Vl^2 
97(2-7)-|-4Vl^2 



(2.26) 



In the Brans-Dicke frame these are the power law solu- 
tions for a barotropic fluid in Brans-Dicke gravity as flrst 
given by Nariai |2^j2l[], where the power law exponents 
of the scale factor and scalar field are given by 



Pi 



qi 



2w(2 - 7) -I- 2 
3cj7(2 - 7) -t- 4 

2(4-37) 
3^7(2-7) +4 



(2.27) 



2. Kinetic dominated solutions 



X2,3 = ±1 



y2,3 = 



(2.28) 



These solutions correspond to negligible fluid density and 
negligible scalar field potential so that the Friedmann 
constraint equation is dominated by the kinetic energy of 
the scalar field with a stiff equation of state 7^ = 2. We 



have power-law solutions of the form given in Eq. (2.1^ 
with 



P2,3 = 1/3 



92,3 =±V2/3 



(2.29) 



In the Brans-Dicke frame we recover the vacuum so- 
lutions of Brans-Dicke gravity found by O'Hanlon and 
Tupper [^ , ^ with power-law exponents given by 



P2,3 



92,3 



%/6t3v/2(3 + 2w) 

±2V6 
3^2(3 2w) T \/6 



(2.30) 



3. Scalar field dominated solutions 

XA - A/V6 , ?/4 = (1 - AV6)1/2_ (2.31) 

Here the fiuid density is negligible, but neither the kinetic 
energy, nor the potential dominates the energy density of 
the scalar field in the Einstein frame. The scalar field has 
an effective equation of state 7^ = A2/3, and for A^ < 2 
this solution corresponds to power-law inflation 0. For 
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A ^ we have power-law solutions of the form given in 
Eq. (|2.18[) with 



P4 = 



A2 



94 = 



A 



(2.32) 



In the Brans-Dicke frame the power law exponents are 
given (for n ^ 1 and n 7^ 2) by ||27|| , 



Pi 



2uj + n+l 
(2 - n)(l - n) ' 
2 



94 = 



1-n 



(2.33) 



For n — the potential acts like a false vacuum energy 
density in the Brans-Dicke frame and this corresponds to 
extended inflation solutions for w > 1/2 For n = 2 
it is the potential in the Einstein frame that remains 
constant, leading to de Sitter expansion |2^,Q. The case 
n = 1 was studied by Kolitch |34|. 



4- 3- Way, Matter-Kinetic-Potential scaling solutions 
37 



Xc, = 



VQ\-2W' 
97(2 - 7) - 2M^(V6A - 2W) 



{V&\-2Wf 



(2.34) 



Here neither the fluid nor the scalar field dominates the 
evolution, and we have self-similar solution where both 
the potential and kinetic energy density of the scalar field 
remains proportional to that of the barotropic matter. 
The effective equation of state for the scalar field is given 

by, 



7i>5 = 7 



67 



67 - T^(%/6A - 2W) 



(2.35) 



We have power-law solutions of the form given in 
Eq. ( ^l8|) with (l|Jl| 



2_ ( V6A-2P^ 
3^ \ T^A 



95 



(2.36) 



In the Brans-Dicke frame this corresponds to power- 
law evolution with exponents given by. 



37 \n — 1 



95 



1 — n 



(2.37) 



As far as we are aware, this solution has not been dis- 
cussed before in the context of Brans-Dicke gravity. It 
is interesting to note that the cosmological evolution in 
the Brans-Dicke frame is independent of the Brans-Dicke 
parameter although it does determine the existence of 
this 3- way scaling solution (see Appendix). 



C. Recovering the Einstein limit when uj ^ 00 

There are two limiting procedures which one can ap- 
ply to recover minimally coupled general relativistic so- 
lutions. 

One can work with quantities defined with respect to 
the Einstein frame where we always have a fixed gravita- 
tional coupling Gm- We can take the limit — ^ to de- 
couple the barotropi c flu id fro m th e scalar fleld but keep 
A flnite. From Eqs. (2.7) and (2.17) this requires lo ^ 00 
but also |n| ^ 00 in the Brans-Dicke frame. In this 
way we obtain a minimally coupled scalar field in Ein- 
stein gravity with an exponential self-interaction poten- 
tial which possesses self-similar solutions where the en- 
ergy density of the field scales with that of the barotropic 
matter for A^ > 37 [^. However in this case the scalar 
field never changes the cosmological dynamics from that 
in a conventional universe dominated by the barotropic 
fluid on its own, so these solutions by themselves can 
never produce an accelerating universe at late times. 

Alternatively one can let w — > 00 in the Brans-Dicke 
frame while keeping G = ijjj(2-K(\?\ ^^(0) and n fi- 
nite. This means that finite changes in the scalar field 
(/) ^ (/) -|- leave both the gravitational coupling G and 
the potential energy V fixed. This corresponds to the 
familiar general relativistic solution for a massless scalar 
field plus cosmological constant. The cosmological con- 
stant will always dominate over any barotropic fiuid with 
7 > 0, but there is always a fine-tuning required if we are 
to postpone this domination until the present epoch. 

In the context of a scalar-tensor gravity theory it is 
intriguing to note the implicit connection between fixing 
the gravitational coupling and fixing the effective value 
of the cosmological constant. A theory in which uj has 
small effective value in the early universe can naturally 
give rise to a self-similar cosmological evolution where 
the potential energy scales with radiation or dust. If lo 
itself can dynamically evolve to large values and fix the 
gravitational coupling at late times [ p5[ then it would also 
fix the value of the effective cosmological constant. But 
a variable ti-'(0) theory breaks the scale-invariance of the 
pure Brans-Dicke theory and requiring that becomes 
fixed only at late-cosmological times is another manifes- 
tation of the cosmological coincidence problem [|| . 



III. COSMOLOGICAL QUINTESSENCE 

We will now consider whether the self-similar solutions 
presented above could form the basis of a viable cos- 
mological model which produces an effective cosmolog- 
ical constant dominating only at late-times without fine- 
tuning. There are two key features which appear promis- 
ing in this respect. Firstly the scalar field is automat- 
ically decoupled, = 0, for relativistic matter (radia- 
tion) with a traceless energy-momentum tensor, 7 — 4/3. 
In this case the Einstein frame evolution coincides with 
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that for a minimally coupled scalar field with exponen- 
tial potential which has scaling solutions for > 87. 
Secondly, we naturally expect a change in the cosmolog- 
ical evolution when non-relativistic (pressureless) matter 
with 7 = 1 comes to dominate the density at late times. 
Here the non- minimal coupling could play a key role and, 
at least in principle, it seems possible to effect a change 
to an accelerating universe. 

These two phases of evolution are fundamentally dif- 
ferent and the change from radiation to matter domi- 
nation ("matter-radiation equality") is known to occur 
relatively late in the cosmological history. We consider 
separately these two phases of cosmological evolution. 

In a multi-fluid cosmology we need to consider the pos- 
sibility that different fluids are minimally coupled with 
respect to different metrics. Relativistic matter (radi- 
ation) is minimally coupled in any conformally related 
frame, but in the case of non-relativistic matter we will 
distinguish two alternative scenarios: 

1. Brans- Dicke theory: In the original gravity theory 
proposed by Brans and Dicke it assumed that 
all matter was minimally coupled with respect to 
the Brans-Dicke frame metric g™. This respects 
the weak equivalence principle |3q |. Thus all mat- 
ter, ordinary baryonic matter and any dark matter, 
must be minimally coupled to scalar field, in 
the Brans-Dicke frame. Conformally transforming 
to the Einstein frame introduces an explicit energy 
transfer to the scalar field, ip, [given by dQ/dt in 



Eq. (2.12)] for all non-relativistic matter. 



2. Generalised dark matter: Given the uncertain form 
of the cold dark matter assumed to dominate the 
matter density in the universe today, we will also 
consider the possibility that this dark matter and 
ordinary baryonic matter are minimally coupled 
with respect to different metrics p7| , ^8|Jl2| . This 
violates the weak equivalence principle but experi- 
mental tests are only performed with baryonic 
matter and baryonic observers. We will consider in 
what follows the case where baryonic matter is min- 
imally coupled in the Einstein frame {dQ-e/di — 0). 
This is equivalent to setting wb — *■ 00 as far as bary- 
onic observers are concerned, which is consistent 
with weak-field tests But some or all of the 

dark matter is minimally coupled in a Brans-Dicke 
frame with finite ojdm which leads to dQ^M/dt 7^ 
in the Einstein frame. 

Amendola considered case (1) where observational 
constraints on the value of ui restrict the cosmological 
evolution to be close to the conventional general relativis- 
tic case. Instead we will consider the second possibility 
of generalised dark matter where solar system tests of 
weak-field gravity do not directly constrain wdm- It is 
important that in this case cosmological observations of 
baryonic matter refer to the quantities calculated in the 
Einstein frame. 



A. Radiation 

The late-time attractor for a decoupled barotropic fiuid 
{W = 0) plus scalar field with an exponential potential 
in the Einstein frame is either the scalar field dominated 



solution given in Eq. (|2.32|) for flat potentials with < 



87 or the 3-way scaling solution given in Eq. ( ^.36 ) for 



steeper potentials with A > 87. This latter solution 
yields acuP where p = 2/(87) ^ previously studied [^,1). 

The energy-momentum tensor for radiation with 7 — 
4/8 has vanishing trace, so it is automatically decou- 
pled from a scalar fleld non-minimally coupled to the 
spacetime scalar curvature, as shown by Eq. (2.17). The 
only constraint that can be placed upon the evolution of 
the fleld during this radiation era is that the scalar field 
should not disrupt the standard temperature-time rela- 
tion at nucleosynthesis. If ordinary barotropic matter, 
specifically protons and neutrons, as well as electrons and 
neutrinos are minimally coupled in the Einstein frame, 
then this constraint must be applied in this frame where 
p = 1/2 is unperturbed, and we simply require that the 
scalar field does not contribute a significant fraction of 
the overall energy-density. Substit uting Eq. (2.84) for 
the 8-way scaling solution into Eq. ( 2.21 ) we have 



A^ > ^ 



^Max 



(8.1) 



where the current bound on the energy density of the 
scalar field at the time of nucleosynthesis, ^^'^^.^ is in 
the range 0.08 to 0.2 |3|,|§. 

For large |A| the steepness of the effective potential, 
ensures, rather paradoxically, that the evolution of the 
scalar field tj) is slow. This is because the gradient 
of an exponential potential depends exponentially upon 
the value of the scalar field. As the scaling solution 
is approached the gradient of the scalar field is bal- 
ance d aga inst the frictional effect of the radiation density 
Eq. (2.10). The larger the exponent, the slower the scalar 
field must evolve in order to maintain scaling. From 
Eq. ( |2.7D we see that the limit |A| ^ 00 can be obtained 
either as |n| — > 00 for lo > —3/2 or for finite n when 
Lu — > —8/2. This is in marked contrast to the more fa- 
miliar case of nucleosynthesis constraints on Brans-Dicke 
gravity where we require weak-coupling lu ^ 00 ml 42 1 . 



B. Non-relativistic matter 

Most estimates of the present density of non- 
relativistic (pressureless) matter in the universe are in 
the range O^.o = 0.2 - 0.4 @|4|. There are upper fim- 
its on the fraction of baryonic matter coming from mod- 
els of primordial nucleosynthesis. Observational upper 
bounds on the abundance of deuterium in the universe 
places an upper limit on the baryonic matter content, 
51b/i^ — 0.02 En], leaving most of the matter as some 
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unknown cold dark matter. We consider here the cos- 
mological evolution where this CDM is non-minimally 
coupled to the scalar field ijj in the Einstein frame. 

In addition recent observations imply that at the 
present time the universe is undergoing an accelerated 
expansion 46 1, which equates to p > 1 in the Einstein 
frame for power-law solutions presented in Sect. II B. 



This is only possible for the scalar field dominated so- 



lution given in Eq. (2.32) if A < 2, which would lead to 



acceleration even during the radiation era which is wildly 
incompatible with primordial nucleosynthesis, or for the 
3-way matter-kinetic-potential scaling solution given in 
Eq. ( 2.36| ) for A < 0. In practice the nucleosynthesis con- 
straint on A^ during the radiation era, given in Eq. , 
ensures that the 3-way scaling solution must exist and is 
the unique late-time attractor for non-relativistic matter 
(7 = 1) and A < 0. 

Acceleration, ps > 1, for the 3-way matter-kinetic- 
potential scaling solution with 7 = 1 and A < requires 



\/6A < AW 



(3.2) 



which from Eqs. ( |2.7|) and ( 2.17 ) implies tha t 2 < n < 
3 in the Brans-Dicke field potential in Eq. ( |2.3| ). The 
existence condition (see Appendix) reduces to 



V6A > -W + y/W^ + 18 



which, combined with Eq. (3.2) gives as a necessary 
condition for an accelerated 3-wa y scaling solution that 
W > v/3/2 which, from Eq. (|2.17|) , impHes w < 1/2. Be- 
cause nucleosynthesis places such a strong lower limit on 
the allowed value of A^ this requires the coupling con- 
stant Lo, to be close to its minimum value of —3/2 for 
2 < n < 3. 

The present day density para meter for t he scalar field 
-0 is then given from Eqs. ( 2.2l| ) and (2.34) by. 



2W{2W ~ V6A) 18 
{2W - V6A)2 

2n + \h + 12^DM 
(2n-3)2 ■ 



(3.4) 



Requiring an accelerating universe at the present day 
places a lowe r bound on the coupling W for a given A 
[see Eq. (3.2)] and satisfying the nucleosynthesis bound 
on A during the radiation era places a lower bound on 
the value of the present day energy density of the scalar 
field V': 



1 



■0,0 



> 



^Max 

ip^nuc 



(3.5) 



Thus a small energy density in the scalar field during the 
radiation era, consistent with nucleosynthesis bounds, is 
compatible with a large density in the present due to 
the different couplings of relativistic and non-rclativistic 
matter to the scalar field. 




0.8 0.7 0.6 0.5 0.4 0.3 0.2 

Present day matter density, Hfyi q 

FIG. 1. QMfi-Wq parameter space for 3-way mat- 
ter-kinetic-potential scaling solution with 7 = 1 and 
A < 0. Limits from nucleosynthesis constraints requiring 



t/j ,nuc 



— 0.2 are shown by a dot-dash line. The solid line 
shows the constraint from requiring acceleration at present. 

The scalar field behaves like a decaying cosmological 
constant whose effective equation of state can be de- 
scrib ed by Wq = — 1, which is given from Eqs. (2.22) 
and (|2.34| ) as 

(3.e) 



Wa 



2W - V6A 2n - 3 



(3.3) The relationship between Cl 



M,0 



1 



■0,0 



and 



Wn IS 



shown in Figure 1. In the limit of strong coupling be- 
tween the scalar field and the dark matter {W ^ 00 or 
i-^DM ^ ^3/2) the equation of state approaches that of a 
cosmological constant — > — 1. In contrast to the usual 
cosmological models where a late time effective cosmolog- 
ical constant is recovered only for a sufficiently flat poten- 
tial, here |A| — *■ 00 is allowed, whilst in the conformally 
related Brans-Dicke frame n remains finite [see Eq. (2.7)]. 
We see from Figure 1 that if we require flM,o ~ 0-3, then 
the nucleosynthesis constraint on |A| drives Wq very close 
to -1. 



IV. CONCLUSIONS 

We have presented the results of a phase-plane analysis 
of the evolution of spatially flat FRW cosmological mod- 
els containing a barotropic fluid in a scalar-tensor theory 
of gravity where the Brans-Dicke scalar, oc 0^, has a self- 
interaction potential, V oc 0^". In contrast to minimally 
coupled flelds in general relativity, this scalar held has 
self-similar scaling solutions for an arbitrary power-law 
potential which correspond to flxed points in the phase- 
plane. The scalar-tensor gravity theory with perfect fluid 
in the Brans-Dicke frame, is equivalent to general rela- 
tivistic gravity with an explicit interaction between the 
fluid and the scalar field in the conformally related Ein- 
stein frame. In this frame the effective potential for the 
scalar field is an exponential and the phase-plane coin- 
cides with that studied recently by Amendola |l4|, and 
Bfllyard and Coley |l|]. 
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There are five different possible fixed points for the sys- 
tem. Two of these are solutions where the scalar field's 
kinetic energy dominates, and the fixed points coincide 
with O'Hanlon and Tupper's vacuum solutions in Brans- 
Dicke theory A third solution, where the scalar 

field's kinetic and potential energy remain proportional, 
coincides with Mathiazhagan and Johri's |Q scalar-field 
dominated power-law solution, which is inflationary for 
< 2. The fluid density is non-negligible only at the two 
remaining fixed points. When the self-interaction poten- 
tial energy of the scalar field is negligible, we recover Nar- 
iai's power-law solutions for a barotropic fiuid in Brans- 
Dicke gravity. In these solutions, the rate of change of 
the Brans-Dicke field, and hence the effective gravita- 
tional constant G oa 1/0^, is inversely proportional to 
the Brans-Dicke parameter, lo. But there is also a new 
solution when the scalar field kinetic and potential energy 
densities both remain proportional to the fluid density. 
In contrast to Nariai's solutions, the rate of change of 
the effective gravitational constant is independent of the 
value of u and depends only upon the potential exponent 
n. This 3-way scaling solution was originally discovered 
by in the context of Einstein gravity ||lj,^ . 

We have investigated whether it is possible in this the- 
oretical framework to construct cosmological solutions 
with non-relativistic (pressureless) matter which are ac- 
celerating at late cosmic times. A generalised dark mat- 
ter model may be one way to achieve this, where non- 
relativistic dark matter is minimally coupled in a Brans- 
Dicke frame, but baryonic matter is minimally coupled in 
the conformally related Einstein frame (consistent with 
the current bounds on the variability of the gravitational 
constant |4^). This leads to an explicit energy transfer 
between the dark matter and the scalar field in the Ein- 
stein frame. In contrast to other scalar field models of 
quintessence [^p^] , scalar fields with potentials of pure 
exponential form, but non-minimally coupled to the dark 
matter, allow a late-time scaling solution with equation 
of state Wq ~ —1, as shown in Figure 1. It is only possi- 
ble for the scalar field to make a significant contribution 
to the present energy density, while remaining consistent 
with nucleosynthesis limits upon the evolution during the 
radiation era, by approaching the strong coupling limit 
(wdm — ~3/2). Such strongly coupled dark matter is 
very different from the usual model of non-interacting 
dark matter and a more detailed investigation is neces- 
sary to determine whether this could represent a viable 
cosmological model 0. 

Due to the scale-invariant nature of the theory and the 
existence of self-similar cosmological solutions, there is no 
fine-timing required of the scalar field, either its initial 
conditions or its parameters, in order to produce a change 



■'■Amendola has recently investigated large-scale struc- 
ture constraints on a coupled quintessence model. 



to an accelerated expansion at late cosmic times. This 
can occur due to a change in the cosmological equation 
of state when a non-relativistic dark matter component, 
directly coupled to the scalar field, begins to make a sig- 
nificant contribution to the total density. This is a fa- 
miliar feature of the standard hot big bang model where 
non-relativistic species, whether baryons or massive neu- 
trinos or more exotic dark matter candidates, only come 
to affect the expansion at relatively late cosmic times. 
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APPENDIX: APPENDIX: EIGENVALUES OF 
CRITICAL POINTS 

To determine the stability of solutions we expand 
about the critical points {xi,yi), 



X % I 



(A.i) 



which to first order gives the equations of motion for u 
and u, 



0- 




((AV6 


- ilXi - 


(3(2- 






- - 



We solve the simultaneous first order equations for two 
eigenvalues m± such that the eigenmodes u — K±v 



evolve as 



Stability of the critical points 



(xi, yi) depends on the eigenvalues m±. Points are stable 
if the real part of the eigenvalues are negative, linear per- 
turbations about the point being exponentially damped. 
A saddle point in the phase plane has the real part of the 
eigenvalues of opposite sign, whilst instability requires 
both real parts to be positive. The eigenvalues for the 
system are listed below. 

In what follows we will make extensive use of the fol- 
lowing bifurcation values for the parameters: 



(A.3) 
(A.4) 
(A.5) 
(A.6) 



7a 




5 




= \/6 , 






97(2- 


7) 


VGXd± 


= w± ^/w^ + 


187 
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1. 2-Way, Matter-Kinetic scaling solution 



[4] I. Zlatev, L. Wang and P. J. Steinh ardt, Phys. Rev 



TO- 



(2-7)2-(2-7a)= 



m+ = — - 



(2-7) 
2W{X-Xc) 
V6(2-7) 



(A.7) 



Exists for 7 < 7a and stable for W{X — A^) > 0. 
2. Kinetic dominated solutions a; = ±1,?; = 



TO+ ==3(2-7)=f21¥, 



(A, 



always exist, x — ±1 stable for ±W^ > 0, 7 > 7a 
and ±A > A^. 

3. Scalar field dominated solutions 



to+ = --(6-A2) , 

m_ = (A - A<j-)(A - Xd+) 



(A.9) 



Exists for A^ < A^ and stable for Xd- < A < Xd+- 

4. 3- Way, Matter-Kinetic-Potential scaling solutions 

Eigenvalues, m± are defined by the solution to the 
quadratic equation. 



+ bm + c = 0. 

where b and c are defined as, 
3 



(A.IO) 



b = 



{V6x~2wy 



[3(2-7)(A-Ad_)(A-Ad+) 
-2WV6{X - Ac)] , 



-6V6W{X - Xd-)iX - Xd+)iX - Ac) 
(V6A - 2WY 



(A.ll) 



Exists for Vl^(Ac — A) > and either A < A^- or 
A > Xd+- Always stable when it exists. 
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